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ABSTRACT 

We estimate cluster masses and velocity dispersions for 123 clusters from optical spectroscopy to 
compare the Sunyaev-Zeldovich (SZ) mass proxy and dynamical masses. Our new survey, HeCS-SZ 
(Hectospec Cluster Survey of SZ-selected clusters), includes 7,721 new or remeasured redshifts from 
MMT/Hectospec observations of 24 SZ-selected clusters at redshifts 2:=0.05-0.20 and not in previous 
surveys. We supplement the Hectospec data with spectra from the Sloan Digital Sky Survey (SDSS) 
and cluster data from the Cluster Infall Regions in SDSS (CIRS) project and the Hectospec Cluster 
Survey (HeCS), our Hectospec survey of clusters selected by X-ray flux. We measure the scaling 
relation between velocity dispersion and SZ mass estimates from the integrated Compton parameter 
for an SZ complete sample of 83 clusters. The observed relation agrees very well with a simple virial 
scaling from mass (based on SZ) to velocity dispersion. The SZ mass estimates (calibrated with 
hydrostatic X-ray mass estimates) are not significantly biased. Further, the velocity dispersion of 
cluster galaxies is consistent with the expected velocity dispersion of dark matter particles, indicating 
that galaxies are good dynamical tracers (i.e., velocity bias is small). Significant mass bias in SZ 
mass estimates could relieve tension between cosmological results from Planck SZ cluster counts and 
Planck CMB data. However, the excellent agreement between our measured velocity dispersions and 
those predicted from a virial scaling relation suggests that any SZ mass bias is too small to reconcile 
SZ and CMB results. In principle, SZ mass bias and velocity bias of galaxies could conspire to yield 
good agreement, but the required velocity bias is agaiaxy ~ 0.77cr£)M, outside the range of plausible 
models of velocity bias in the literature. 

Subject headings: galaxies: clusters: individual — galaxies: kinematics and dynamics — cosmology: 
observations 


1. INTRODUCTION 

As the universe evolves, the comoving number den¬ 
sity of clusters of fixed mass increases. The evolu¬ 
tion of cluster abundances depends strongly on the 
amount of dark matter and dark energy in the uni¬ 
verse. Thus, many groups have used different clus¬ 
ter mass proxies to determine the mass f unction and 
const r ain cosmological parameters fe.g.. iR ines et ^ 
2f)f)7l IWia IVikhliriin et^ l2nn9H: IHenrv et all I2nnft 
Mantz et al.l I2010IJ : iRozo et al.1 120101 and references 
therein). Recently, o thers have used the Suny aev- 
Zeldovich (SZ) effect (|Sunvaev fc Zeldovichl 1197^ to 


identify large sam ples of clusters to constrain cosmo- 


logical parameters (IBenson et al. 


Hasselfield et al.1 

120131 IPlanck Collaboration et al. 

2014b 



Data from the Planek satellite show that cosmologi¬ 
cal parameters determined from anisotropies in the cos¬ 
mic microwave background disagree with those derived 
from cluster abunda nce measurements from the Planek 
SZ cl uster survey ((Planck Collaboration et ahl l2014aL 
I2015a() . Fewer clusters are observed than predicted by 
the cosmology that best fits the Planck CMB data. Inter- 
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estingly, estimates of the amplitude of structure from cos- 
mic shear yield a simila r tension with Planck CMB data 
((MacCrann et al.l 1201511 . If SZ masses (calibrated from 
X-ray observations) systematically underestimate true 
masses by about 45%, the cosmological parameters de¬ 
rived from SZ cluster counts shift into agreement with the 
CMB results (jPlanck Collaboration et al.ll2014;ih . An al¬ 
ternate analysis using weak lensin g data for mass cali¬ 
brati o n finds no significan t tension ( von der Linden et al.l 
12014 iMantz et al.l 1201511 , suggesting that the tension 
could arise from biases in the calibration of SZ masses. 

Here, we compare SZ mass estimates to dynamical 
mass estimates based on the redshifts of cluster members. 
Dyna mical m a ss est i mates have a long history beginning 
with iZwickvl (119331 . 1193711 . In numerical simulations, 
either the virial theorem or the caustic technique can 
provide cluster mass estimates with little bias but with 


some 

intrinsic scatter dne to projection effects (IDiaferio 

1999 

Evrard et al.ll2008HSerra et al.ll2011tlMamon et al. 

2013 

Gifford & Milleill2013HOld et al.ll2014ll. Hvdrodv- 


namical simulations show that the velocity distribution 
of gal axies is very similar to that of dark matter par - 
ticles ((Faltenbaclier fc DiemancU 120061: iLau et ahl 1201011 , 
with the possible exception of the brightest few galaxies 
(|Lau et al.1 120101: IWu et al.1 1201311 . Thus, virial masses, 
caustic masses, or dynamical mass proxies such as veloc- 
it y dispersion a.re a p owerful test of SZ mass estimates. 

iRines et al.l (|2010ll made the first comparison of SZ 
signals to mass estimates from galaxy dynamics, but 
the sample was limi ted to 15 clusters. A later study 
bv ISifon et al.1 (1201311 obtained optical spectroscopy for 
16 SZ-selected clusters selected from observations with 
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the Atacama Cosmology Telescope (ACT); they found 
that the scaling relation between SZ signal and mass 
(actually measured from velocity dispersions) is consis¬ 
tent with relations d etermined with ot her mass calibra¬ 
tors (X-ray, leasing). iR.uel et id] (j2014ll measured veloc¬ 
ity dispersions for SZ-selected clusters identified in ob¬ 
servations with the South Pole Telescope (SPT); they 
conclude that SZ signal co rrelates well with vel ocity dis¬ 
persion. The SPT results (iBocauet et al.l[2M^ are con¬ 
sistent with positive velocity bias (that is, the velocity 
dispersion of the galaxies is larger than the velocity dis¬ 
persion of the dark matter particles). The clusters in the 
ACT and SPT samples span a wide range of redshifts 
(0.2 < z < 1.3). It is possible that the scaling between 
velocity dispersion and virial mass evolves significantly 
over that period. Further, the spectroscopy for these 
clusters is often incomplete at large radii or contains rel¬ 
atively few cluster members. In principle, the measured 
veloc i ty dispersions c ould be biased fe.g.. iBiviano et ahl 
l2nMlWu etHllml . 

To provide a much broader foundation for compari¬ 
son of dynamical and SZ mass proxies, we compare SZ 
mass estimates of 123 clusters from the Planck SZ cata¬ 
log with velocity dispersions from wide-held optical spec¬ 
troscopy. Several clusters have reds hifts in the Sloan Dig¬ 
ital Sky Survey (lAhn et al.l l2014fl , and many are part 
of the Cluster Infall R egions in SDSS project (CIRS; 
IRines fc Diaferioll2006n or the Hectospec Cluster Survey 
iHeCS: IRines et al. [1^. To supplement this sample 
and create an SZ-selected sample of clusters, we con¬ 
ducted HeCS-SZ, an MMT/Hectospec spectroscopic sur¬ 
vey of 24 clusters. We also include analysis of 30 clusters 
from SDSS redshifts. 

We discuss the cluster samples and spectroscopic data 
in §2. We measure the SZ-optical scaling relations in 
§3. We discuss the implications of our results in the 
context of other cosmological observations in §4. We 
assume a cosmology of Dm=0.3, Da=0.7, and Ho=70 
km s“^ Mpc“^ for all calculations. 

2. OBSERVATIONS 

2.1. Optical Photometry and Spectroscopy 

HeCS-SZ is an extension of the HeCS survey to include 
clusters that enable construction of an SZ-limited sam¬ 
ple. We observed 7,721 new redshifts in 24 clusters. We 
combine these new measurements with the existing HeCS 
and CIRS surveys and with data from the literature to 
construct a total sample of 123 clusters. For all but a few 
clusters the sampling is sufficient for a robust determi¬ 
nation of velocity dispersion. We use SDSS photometry 
for all clusters. 

2.1.1. Spectroscopy: CIRS and HeCS 

The Hectospec Cluster Survey (HeCS) is a spectro¬ 
scopic survey of 58 galaxy clusters at moderate redshift 
(z=0.1-0.3) with MMT/Hectospec. HeCS includes all 
clusters with ROSAT X-ray fluxes of /x > 5 x 10“^^erg 
s~^ at fO.5-2.0] keV from the Bright Clus ter Survey (BCS; 
Ebeling et al.lll99^ or REFLEX survey (iBdhringer et alJ 
2004f) with optic al imaging in the Sixth Data R elease 
(DR6) of SDSS (jAdelrnan-McCarthy et al.ll2008[) . We 
used DR6 photometry to select Hectospec targets. The 
HeCS targets are all brighter than r=20.8 (SDSS cata- 



Fig. 1. — Planck SZ mass estimates versus redshift. Open black 
squares show clusters with dynamical mass estimates from CIRS. 
Solid blue and red squares show clusters from HeCS and new HeCS- 
SZ clusters respectively. Small points show the remainder of the 
Planck SZ catalog. The clusters studied here are representative of 
clusters at z<0.3 in the Planck SZ catalog. 

logs are 95% complete for point sources to rfv22.2). 

For He CS, we acquired spectr a with the Hectospec in¬ 
strument (|Fabricant et al.l[2005r) on the MMT 6.5m tele¬ 
scope. Hectospec provides simultaneous spectroscopy of 
up to 300 objects across a diameter of 1°. This telescope 
and instrument combination is ideal for studying the 
virial regions and outskirts of clusters at these redshifts. 
Because cluster properties such as projected velocity dis¬ 
persion depend on radius, wide-held spectroscopic cover¬ 
age is important for measurin g accurate global ve locity 
dispersions and virial masses (IBiviano et al.ll200^ . We 
used the red sequence to preselect likely cluster members 
as primary targets, and we Filed otherwi se unassigned 
hbers with bluer targets (|Rines et al.ll2013L describes the 
details of target selection). 

CIRS used spectroscopy from the Fourth Data Release 
of SDSS to study the virial and infall regions of clusters. 
We use the dynamical data tabulated in CIRS for 25 
clusters. We update dynamical parameters for two addi¬ 
tional CIRS clusters: A2249 was poorly sampled in DR4 
but has many more redshifts available in DRIO. We use 
the DRIO redshifts to update the dynamical parameters. 
The central region of A2175 was poorly sampled in DR4. 
We thus obtained additional redshifts in the central parts 
of A2175 with Hectospec (see below). 

2.1.2. Spectroscopy: HeCS-SZ 

We observed 24 clusters in the Planck catalog of SZ 
clusters using MMT/Hectospec. The target clusters are 
in the redshift range 0.05 < z < 0.20 and were observed 
mostly in decreasing Planck signal-to-noise ratio (a few 
clusters with relatively weak SZ signals were observed as 
backup targets for variable observing conditions). We 
also observed one held in A2175, a cluster from CIRS 
with limited SDSS spectroscopy in CIRS. Preliminary 
analysis indicated that A2175 had an unusually small 
velocity dispersion given its SZ mass. The additional 
redshifts in A2175 show that the CIRS data led to a 
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Fig. 2. — ROSAT 'K.-ray luminosities of Plaracfc-selected clusters 
versus redshift. Filled symbols are clusters in the HeCS-SZ sample: 
black squares show clusters with dynamical mass estimates from 
CIRS, blue points are clusters from HeCS, and red points are new 
clusters in HeCS-SZ. 

significant underestimate of its velocity dispersion and 
caustic mass. 

Our observing strategy closely matches that of HeCS: 
we used SDSS photometry to identify a red sequence in 
each cluster field. We then identify a cutoff in appar¬ 
ent magnitude that offers a good compromise of high 
completeness (sparser targets produce fewer fiber con¬ 
flicts) and dense sampling. Targets are primarily drawn 
from galaxies with g — r colors within 0.2 mag of the 
red sequence, and we assign higher priorities to brighter 
galaxies and galaxies closer to the cluster center. This 
approach provides reasonably high sampling in the clus¬ 
ter cores but can lead to relatively sparse sampling of 
dense regions outside the core. We included galaxies with 
slightly bluer colors (up to 0.4 mag bluer than the red 
sequence) as targets to fill fibers when available. We 
matched all targets to redshifts from the literature as 
compiled by NEE0 as of 2013 September as well as to 
SDSS DR8 spectra. Most of the targets with existing 
redshifts are from SDSS, but seve ral are from targeted 
studies of individual clusters ('e.g.. iCvoriano et "311120051 
for A586). Targets with existing redshifts are removed 
from the targeting catalogs prior to fiber assignment. 

Table [T] lists 7,721 new redshifts measured with Hec¬ 
tospec. We visually inspected all spectra to confirm the 
reliability of the redshift. Column 5 of Table [T] lists 
the cross-c orrelation score Rxn from the IRAF pack¬ 
age rvsao (|Kurtz fc MinkI [r998ll . A score of Rxc > 3 
indicates a reliable redshift; some galaxies with smaller 
values of Rxc are included when visual inspection shows 
multiple obvious absorption and/or emission lines and 
the spectrum suffers from contamination (e.g., light 
bleeding into the spectrum from a nearby fiber contain¬ 
ing a bright star). Table [2] lists redshifts from SDSS and 
other literature (as compiled by NED) for galaxies classi¬ 
fied as cluster members by the caustic technique (see be¬ 
low). Table|3]lists 168 redshifts measured with the FAST 

® http://ned.ipac.caltech.edu 


instrument (|Fabricant et al.lll998ll on the 1.5-meter Till- 
inghast telescope at the Fred Lawrence Whipple Obser¬ 
vatory. The additional single-slit spectra from FAST re¬ 
duce the incompleteness of bright (SDSS r < 16.5) galax¬ 
ies in the HeCS-SZ clusters. 

In addition, we identified several clusters in the Planck 
SZ catalog that lie below the completeness limits but that 
are at sufficiently low redshift {z < 0.1) that they have 
reasonable redshift coverage in SDSS DRIO. We include 
these clusters in an extended sample. 

We include four nearby {z < 0.05) clusters that lie 
inside the SDSS DRIO photometric footprint but outside 
the SDSS spectroscopic footprint. These nearby clusters 
have large numbers of redshifts available in the literature. 
Because of the redshift dependence of the limiting mass 
for SZ detection by Planck (driven by the large beam size 
of Planck), including these low-redshift clusters improves 
the sampling of low-mass clusters in the sample. The 
FAST redshifts in Table [3] are especially useful for these 
clusters. 

Figure [T] shows the Planck SZ mass estimates versus 
redshift. The minimum mass a cluster must have to be 
detected by Planck increases with redshift because the SZ 
signal of lower-mass clusters at higher redshift is diluted 
by the large beam below the sensitivity of Planck. 

The CIRS and HeCS clusters provide a good sam¬ 
pling of the Msz — z distribution, but this distribution 
is possibly biased due to the underlying X-ray selection 
of CIRS and HeCS. Figure [2] shows the X-ray luminos¬ 
ity of clusters in CIRS, HeCS, and HeCS-SZ as a func¬ 
tion of redshift. The clusters we target with Hectospec 
include clusters that lie above the X-ray flux limits of 
CIRS and HeCS but were not in the appropriate SDSS 
photometric footprint and also clusters that have X-ray 
fluxes below the CIRS/HeCS flux limits. Targeting these 
X-ray-faint clusters enables a test of the impact of X- 
ray selection on the scaling relation parameters based on 
SZ and optical properties. The X-ray luminosities are 
measure d in the ROSAT band but from heterogeneous 
sources (lEbeling et al.lll998HBbhringer et aTll200(i 120051 : 

iPiffaretti et al.ll201lHPlanck Collaboration et al.ll2014bfl . 

A careful study of the X-ray properties of HeCS-SZ clus¬ 
ters would require a homogeneous reanalysis of ROSAT 
X-ray images. 


TABLE 1 

HeCS-SZ Redshifts from MMT/Hectospec 


Coordinates (J2000) 

RA DEC 

C20 

km/s 

<^CZ 

km/s 

Rxc 

Flag 

Member 

00:09:33.60 

32:31:03.16 

83460 

52.19 

6.23 

Q 

0 

00:09:35.55 

32:14:05.00 

69478 

182.17 

1.64 

Q 

0 

00:09:39.32 

32:21:22.38 

122251 

49.89 

5.08 

Q 

0 

00:09:42.73 

32:16:05.45 

83482 

100 

4.45 

Q 

0 

00:09:43.80 

32:33:54.17 

108270 

9.69 

15.29 

Q 

0 


Note. — Table[T]is published in its entirety in the electronic edi- 
tion of the Journal. A portion is shown here for guidance regarding 
its form and content. 


The caustic technique (jDiaferiolll999l[SCTra et al.ll2011[l 
uses a redshift-projected radius diagram to isolate clus¬ 
ter members from foreground and background galax¬ 
ies in phase space. After smoothing the galaxy dis¬ 
tribution in the redshift diagram, the infall regions of 
clusters produce well-defined envelopes containing the 
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Fig. 3. — Redshift (rest-frame clustrocentric velocity) versus pro¬ 
jected radius for galaxies around HeCS-SZ clusters. The caustic 
pattern is evident as the trumpet-shaped regions with high den¬ 
sity. The solid lines indicate our estimate of the location of the 
caustics in each cluster. Clusters are ordered left-to-right and top- 
to-bottom by decreasing mass as estimated from the Planck SZ 
data. 
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Fig. 4.— Same as Figure [3] 


TABLE 2 

HeCS-SZ Members from Literature Redshifts 


Coordinates (J2000) 

RA DEC 

km/s 

O'cz 

km/s 

Ref. 

0:11:45.24 

32:24:56.17 

30309 

100 

2 

0:11:19.72 

32:17:09.39 

32168 

201 

2 

0:20:02.98 

28:44:58.73 

29876 

27 

2 

0:20:05.48 

28:41:01.73 

29545 

47 

2 

0:20:16.85 

28:46:09.69 

26793 

33 

2 


Note. — Table[2]is published in its entirety in the electronic edi¬ 
tion of the Journal. A portion is shown here for guidance regarding 
its form and content. 

Note. — References: [1] SDSS, [2] NED. 


vast majority of cluster members. Specifically, the 


Fig. 5. — Same as Figure[3] 
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Fig. 6.— Same as Figure!^ 


TABLE 3 

HeCS-SZ Redshifts from FLWO 1.5m/EAST 


Coordinates (J2000) 
RA DEC 

C2:0 

km/s 

<^CZ 

km/s 

Rxc 

Member 

0:11:05.08 

31:54:29.53 

24574 

24 

9.01 

0 

0:11:34.79 

32:28:16.28 

30990 

20 

17.46 

1 

0:11:45.24 

32:24:56.20 

30542 

51 

6.47 

1 

0:12:27.58 

32:45:09.84 

12600 

9 

15.40 

0 

0:12:30.47 

32:19:12.45 

24565 

5 

31.67 

0 


Note. — Table[3]is published in its entirety in the electronic edi¬ 
tion of the Journal. A portion is shown here for guidance regarding 
its form and content. 


list of cluster members within r 2 oo Is 96% complete 
and only 2% of the members are actually interlopers; 
within the larger radius 3r2oo , where the caustic tech¬ 
nique is the only usable method, the completeness is 
95% and the interloper fraction is 8% (jSerra fc Diaferiol 
l2013f) . The edges of this distribution are called cans- 
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Fig. 7.— Same as Figure |3] 


tics and they are rela ted to the es c ape velocity pro- 
file of the cluster (see iDiaferiol 120091 : [S erra et ^I 2011 L 
for reviews). The escape velocity profile is the basis 
for a mass profile that can extend into the infall re¬ 
gion where the galaxies are gravitationally bound but 
not virialized. Caustic mass estimates generally agree 
with estimat es from X-ray obse r vations and gravitational 
leasing fe-g.. iRines et alJl2003HBiviano fc Girardil 120031: 


Diaferio et alT200^ Rines fc Diaferiol 120061: iRines et alJ 

2007HGeller et al.ll 201 ,‘lL and references therein). 


Figures [31[7] show the phase space diagrams of the 
HeCS-SZ clusters not already published in CIRS or HeCS 
(the poorly-sampled CIRS clusters A2175 and A2249 are 
reproduced here with enlarged datasets). Almost all clus¬ 
ters display prominent infall patterns, and the caustics 
are shown on the figures. Clusters are ordered by de¬ 
creasing SZ mass, and there is a clear trend of decreasing 
central velocity dispersion with decreasing SZ mass. 

We apply the prescription of iDanese et al.l (|1980ll to 
determine the mean redshift czq and projected velocity 
dispersion Up of each cluster from all galaxies within the 
caustics. We calculate Up using only the cluster members 
projected within r 2 oo estimated from the caustic mass 
profile. Note that our measured velocity dispersions use 
the caustic technique only to define membership and the 
limiting radius r 2 oo- Independent of its performance as 
a mass estimator, the caustic technique is a highly effi¬ 
cient membership selection algorithm, especially at the 
relati vely small radii we focus on here (|Serra fc Diaferiol 
l201.‘lfl . Table m lists the central cluster redshifts, velocity 
dispersions inside r 2 oo, and M 200 from the caustic mass 
profile. The ninth column of Table |4] indicates whether 
the cluster is part of the CIRS, HeCS, or HeCS-SZ sam¬ 
ple. 


2.2. SZ Measurements 

The _ SZ _ measu re ments are from 

iPlanck Collaboration et H] (j2014bf l. an all-sky 
SZ survey. Numerical simulations indicate that 
the integrated Compton y-parameter Ysz has 
smaller scatter t han the peak y-decrement Vnea k 
(|Motl et al.l 120051 : IPlanck Collaboration et ahl l2014bD . 


IPlanck Collaboration et all (l2014bD report only Ysz- 
Although j/peafc should be nearly independent of redshift, 
Ysz depends on the angular size of the cluster. The 
quantity YszD\ remoyes this dependence. Table |4] 
summarizes the Planck SZ measurements. 

The Planck mass estimates are extracted from an aper¬ 
ture of 6 * 500 , the angular radius corresponding to rsoo 
(the radius ta is the radius that encloses a mean density 
of Apc{z) where Pc{z) is the critical density). This ra¬ 
dius is larger than the rad ii probed by some other mass 
estimators. For ins tance. iBonamente et all (|2008f l and 
iMantz et al.l (|2010al l find that X- ray masses are best de- 
termined within r 25 oo 1a lthough iVikhlinin et al.ll2009al 
and others use M^qq). iMarrone et al.l (|2009l l uses an 
aperture of 350 kpc as the best match to their mass esti¬ 
mates from strong grayitational lensing. Because the SZ 
signal falls off more slowly with radius than the X-ray 
flux, the outer parts of clusters are more important for 
SZ obseryables than for X-ray obse ryables. For instance, 
[Planck Collaboration et ahl (l2013bf) used Planck data to 
determine the ayerage pressure profile of the ICM to radii 
of 3 r 5 oo, a regime that is yery difficult to study eyen with 
yery deep Chandra obseryations. Because yirial masses 
and yelocity dispersions are best suited for mass esti¬ 
mates at radii ~ r 2 oo, they may be better suited for 
comparison with SZ mass estimates. 

The central redshifts in the Planck SZ catalog are usu¬ 
ally close to the central redshifts we obtain in our hier¬ 
archical clustering analysis of the cluster redshifts (see 
D99 for details). Howeyer, for about half of the clusters, 
our central redshifts differs by more than a percent from 
the redshifts listed in the Planck SZ catalog. We there¬ 
fore re-scale all SZ integrated Compton parameters by 
[D\{zh)/D\{zsz)] where DA{zh) and Da{zsz) are the 
angular diameter distances for the hierarchical center Zh 
and the Planck catalog redshift zsz- We similarly rescale 
SZ mass estimates using the appropriate scaling relation 
from [Planck Collaboration et al.l (|2014bll . 

We define an SZ-complete sample of 83 clusters from 
the SZ mass proxy M 500 in the Planck SZ catalog. The 
final column of Table U] indicates whether the cluster is 
in the complete sample. The completeness limit corre¬ 
sponds to the 80% completeness limit for the medium- 
deep suryey coyering 44% of the sky and to the 50% 
completeness limit for the s hallow survey coyering the 
remain ing 56% of the sky (jPlanck Collaboration et ^ 
l2014bll . Our sample includes all but four clusters aboye 
this limit: two at moderate redshift (A1677 at 2:=0.18 
and A1759 at z=0.17) and two at low redshift (z w 0.04: 
A2572 and RBS 1929). The SZ completeness limits we 
use are slightly above the 80 % completeness limits of the 
update d Planck SZ catalog (IPlanck Collaboration et ahl 
I2015bll . A quick inspection shows that the updated SZ 
catalog contains few clusters above the completeness lim¬ 
its we use here. 


3. RESULTS 

As discussed in §1, data from the Planck satellite 
indicate tension between cosmological parameters de¬ 
termined from CMB and SZ results. One possible 
resolution to the tension is that the SZ mass esti¬ 
mates (calibrated with hydrostatic X-ray mass esti¬ 
mates) are biased. Comparing dynamical estimates 
of cluster mass from galaxy redshift surveys to the 




















































































Fig. 8.— Scaling relation between projected velocity dispersion <jp and Msz-, the mass proxy based on the integrated Compton parameter 
YsZ' The first two panels show the marginalized probability distribution functions of the parameters a (intercept), b (slope), and Uint (the 
intrinsic scatter). The first panel shows that the slope and intercept are correlated. Shaded regions indicate 68%, 95%, and 99% confidence 
regions from darkest to lightest shading. The crosses show the median and the 68% confidence range of the probability distribution 
functions. The third panel shows the the relation with the median values of a and 5 as a solid line with the median value of the intrinsic 
scatter (Tint shown by the dashed lines. 


TABLE 4 

Dynamical Masses and SZ Signals 


Cluster 

a 

5 

Z 

CTp 

M200,c 

Msz 

YszD\ 

Spectra 

Planck ID 

Sample 


deg 

deg 


km 

IQI^Mq 


10“®Mpc-2 




A0007 

2.93500 

32.41700 

0.10302 

783+“ 

2.77 ± 1.14 

O.Ui/_Q 45g 

U.lUO_Q Q24 

HeCS-SZ 

PSZ1G113.26-29.69 

1 

A0021 

5.17050 

28.67510 

0.09456 


2.92 ± 1.33 

q qoc:~I"0-359 

0.020_q 

U.140_q 025 

HeCS-SZ 

PSZ1G114.78-33.72 

1 

A0076 

10.00200 

6.81800 

0.03999 

4551- 

1.19 ± 0.04 

1-63i12:28 

U.UOZ_Q QQg 

HeCS-SZ 

PSZ1G118.03-55.88 

1 

A0085 

10.45870 

-9.30190 

0.05565 

69212 

2.50 ± 1.19 

4.9001°J2 

n OOK+O-Ois 

CIRS 

PSZ1G115.20-72.07 

1 

A0098S 

11.61470 

20.38645 

0.10380 

59412 

2.17 ± 0.09 

Q 'Tqq-|-0.516 

Z./00 Q 

0 070+04)29 
u.u<y Q Q2g 

HeCS-SZ 

PSZ1G121.35-42.47 

0 


Note. — Table [4| is published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its 
form and content. 

Note. — Redshift 2 and velocity dispersion <7p are computed for galaxies defined as members using the caustics. 


SZ mass proxies tests this hypothesis. Several stud¬ 
ies show a strong correlation betwee n X-ray mass es¬ 
timates and SZ mass estimates (e.g., ^onamente et al. 
2008; lAndersson et al.ll2011I iPlanck Collaboration et ah] 
2011 ; ICzakon et al.l 120151) . but both methods measure 
the properties of the intracluster medium (ICM). Thus, 
systematic effects could still be present. For instance, 
the ICM is likely to depart from hy drostatic equilib¬ 
rium in the outer parts of the cluster (|Bonamente et ^ 
l2013f) . Gravitational lensing does not measure the ICM, 
but it does measure all of the matter along the line of 
sight to the cluster, introdu cing significant scatter into 
lensing mass estim ate s le.g.. iHoekstra et al ll200il 120111 : 

I Hwang et al.l 1201411 . iMarrone et al.l ( 200911 show that 
lensing masses are consistent with SZ estimates although 
with s ignifi cant scatter. R e cently , Ivon der Linden et al.l 
lj2014ll and IHoekstra et all (j2015l l have used large sam¬ 
ples of weak lensing mass estimates to test for systematic 
bias in SZ masses; both groups find that the SZ masses 
are systematically underestimated, but both estimates of 
bias are smaller than the value required to fully recon¬ 
cile Planck CMB a nd SZ results (the uncertain ty range 
in bias obtained bv ivon der Linden et al.l[20l4 includes 
this value within the 2a confidence interval). In contrast, 
iMelin fc Bartlett] (|2014[ 1 use weak lensing of the CMB to 
estimate cluster masses, and they find little evidence for 
mass bias. 

Figure |8] shows the relation between projected veloc¬ 
ity dispersion ap and the mass Msz estimated from the 


Planck data (note that, for most clusters, the measure¬ 
ments of Ysz use X-ray data to determine the region 
where the SZ signal is extracted). We use a Bayesian ap¬ 
proach (see Appendix for details) to determine the best- 
fit relation P{(Tp\Msz), that is, the predicted value of ap 
at a given observed value of Msz- We allow for intrinsic 
scatter in ap that is expected to arise from projection 
effects of non-spherical clusters. Our Bayesian analysis 
yields a relation of 

logioiap) = OM9trolllogio{MYsz) + 2.687t°;°27 (i) 

with cTp in units of km s~^ and Msz in units of 
lO^^h^QMQ. The scatter in ap at fixed Msz is logig ^ = 
0.0973to;oo85- The best-ht parameters and their uncer¬ 
tainties are the medians and the boundaries of the 68% 
conhdence levels derived from the posterior probability 
of the regression parameters. Figure [8] shows this rela¬ 
tion as a solid line. Note that we fit P{ap\Msz) rather 
than the inverse because the statistical uncertainties in 
Msz are smaller than the statistical uncertainties in ap. 

The intrinsic scatter we measure corresponds to about 
a factor of two in the estimated mass within r 2 oo- A 
comparison of several richness-based and dynamics-based 
mass estimators demonstrate similar scatter for several 
mass estimators based on velocity disper sions or vari¬ 
ations of Jeans’ analysis l|01d et al.l 1201311 . Thus, the 
scatter probably represents geometric projection effects 
and not our use of the caustic technique to define cluster 
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Fig. 9. — Scaling relation between projected velocity dispersion 
tTp and the SZ mass proxy Msz based on the integrated Compton 
parameter The thick solid line shows the best-fit relation 

of P{ap\Mgz) with the intrinsic scatter shown as the green band. 
Open squares, filled squares, and open circles represent clusters 
from CIRS, HeCS, and HeCS-SZ r espec tively. The dotted line 
shows the relation from lR,ines et ahl <1201(1 ') from a small number of 
clusters. The red dashed lin e shows the predicted relation using the 
virial scaling relation from lEvrard et al.l <120081 and assuming no 
hydrostatic mass bias (i.e., l-fe=0). The blue dashed line shows the 
same predicted relation for the bias (l-b=0.58) required to reconcile 
Planck SZ counts with the Planck CMB cosmology. 

membership. 

Previous work provides an expected value for this 
slope. Numerical simulations of clusters with a variety 
of codes yield a consistent scaling relatio n of the mass 
M-},nn with velocity dispersion, ap oc M^qq^ (jEvrard et ahl 
1200811 . This slope is measured for randomly selected dark 
matter particles rather than galaxies, but hydrodynami- 
cal simulations suggest that v elocity bias is sin all for large 
samples of cluster galaxies (|Wu et al.l 1201311 like HeCS 
and CIRS (we discuss velocity bias further in il4.4|) . The 
slope of the scaling relation for dark matter particles in 
clusters agrees well with our observed Up — Msz relation 
(Equation 1). Figure |9] shows our data and scaling rela¬ 
tion compared to the virial scaling of dark matter parti¬ 
cles, and the agreement is reasonable. Figure [TUI shows 
the marginalized probability distribution functions of the 
parameters of our scaling relation along with the virial 
scaling of dark matter particles. Figures [S] and [in] also 
show the virial scaling of dark matter particles rescaled 
by assuming that Msz = O.bSMtrue, the mass bias re¬ 
quired to match the SZ counts to the CMB data. Such 
a large mass bias is strongly disfavored by our observa¬ 
tions. 

Figure Hi] shows the best-fit relation for P{ap\YszD\), 
the expect ed velocity dispersion at fixe d SZ ma ss proxy 
YszD\. iPlanck Collaboration et ^ (j2014b[l obtain 
MIqq oc {YszD\) using hydrostatic mass estimates 
from detailed XMM-Newton observations. Because the 
conce ntration-mass relat ion depends weakly on mass 
(e.g., iBullock et al.ll200lll . we use a fixed conversion of 
ATgoo ~ l.SSMsoo appropri ate for concentration c=5 as¬ 
suming an NFW profile ([Navarro et al.l Il997f) . With 
these assumptions, the expected slope of the Cp — Ysz 



b 

Fig. 10.— Parameters of the virial scaling relation between pro¬ 
jected velocity dispersion ap and the SZ mass proxy Msz based on 
the integrated Compton parameter is the intercept, b is 

the slope). Contours show confidence intervals from our Bayesian 
analysis and the cross without a symbol shows the median and 
68% percentiles of the distribution shown with the contour levels. 
Points with errorbars show models based on simulations. The filled 
square is the virial scaling relation of dark matter particles from 
lEvrard et al.l II2Q08I) . The filled circle shows this same relation re¬ 
normalized to reflect a mass bias of Msz = O.bdMtrue, the value 
needed to match SZ and CMB constraints. The other points show 
several models of v elocity bias. The ope n triangle and open square 
show the models of IMunari et al.j 112013) for galaxies identified from 
dark matter subhalos and from hydrodynamical simulations includ¬ 
ing star for mation and AGN feedback. The open diamond shows 
th e model of|Lau et al.l II2010I) , and the open circle shows the model 
of ISaro et al.l ||2013) . 

relation is 0.188 with an intercept of 3.003. Figure [TT] 
shows that these values agree very well with our Bayesian 
analysis (see also Table |S|). __ 

In contrast with our previous work (jRines et al.ll2010ll , 
the relation between projected velocity dispersion Up and 
YszD\ agrees with expectations from scaling relations 
of dark matter particles and simulations of the SZ effect. 
We attribute this difference to both the much larger (5x) 
sample of clusters studied here and the improved statis¬ 
tical methods enabled by the larger sample. 

Figure [T2| shows the best-fit relation P{M 20 [)\Msz), 
the caustic mass M 200 obtained at fixed Msz- The 
intrinsic scatter in this relation is somewhat smaller 
than a factor of two, consistent with the expected scat- 
ter in caust i c ma ss estimates due to projection effects 
(jSerra et al.1 l20lill . Note that a similar level of scat¬ 
ter is foun d for alternate impleme ntations of the caustic 
technique ([Gifford fc Milleill2ni3[ l as well as alternative 
mass estimators based on measured velocity dispersions 
([Old et al.ll2014ll . While a detailed treatment of outliers 
is beyond the scope of this work, we note that one clus¬ 
ter, MS2348-I-2929, with an observed velocity dispersion 
smaller than predicted by its Planck SZ mass, is unde¬ 
tected i n observations with the Arcminute Microkelvin 
Imager ([Perrott et al.ll20I^ . suggesting that the SZ mass 
in the Planck catalog is an overestimate. 

Figure Hni shows the (Jp — Msz relation for the extended 
sample of 123 clusters. There are significantly more out¬ 
liers than in the SZ-complete sample. Most of these out- 
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Fig. 11.— Similar to Figure |8] for the scaling relation between projected velocity dispersion Cp and the integrated Compton parameter 
Ysz- 



Fig. 12. — Similar to Figure |8] for the scaling relation between caustic mass M 200 and Msz, the mass proxy based on the integrated 
Compton parameter Ysz- 


liers have redshifts only from SDSS, and some are at 
z > 0.1. Thus, these clusters are not well sampled. Ob¬ 
taining additional redshifts for these clusters could sig¬ 
nificantly alter the measured velocity dispersions (similar 
to the changes for the CIRS clusters A2175 and A2249 
resulting from additional redshift data from Hectospec 
and SDSS respectively). The best-fit parameters of the 
scaling relation are virtually unchanged, but the inferred 
intrinsic scatter is larger due to the larger number of 
outliers (Table [S|). 

4. DISCUSSION 
4.1. Predictor Relations 

Cluster scaling relations applied to large surveys 
are a basis for cosmological studies, including measur- 
ing the cluster mass function or the power sp e ctrurn 
(|Mantz et al.ll2010E: iRozo et al.ll2010jl . lAndreonl (|2010d 
discusses how, given observable properties A and B, the 
slopes of the predictor relation P{A\B) (the probabil¬ 
ity of a cluster having the property A given an observed 
value of property B) may be significantly different from 
the inverse of the slope of the predictor relation P{B\A). 
This difference is larger when there is significant intrinsic 
scatter in the relation between the two properties. 

Because different investigators require different predic¬ 
tor relations, we include here the relations between sev¬ 
eral mass observables (Tabled. We do not include con¬ 
straints on the Msz — M 200 relation for the extended 
sample because the large scatter caused by a few outliers 
leads to very weak constraints on the parameters of the 



Fig. 13. — Scaling relation between projected velocity dispersion 
(Tp and the SZ mass proxy Msz based on the integrated Compton 
parameter Y^oqD^ for the extended sample of Planck -selected 
clusters (including clusters below the Planck completeness limits). 
Several clusters in the extended sample are outliers below the main 
relation. These clusters are not well sampled in SDSS spectroscopy, 
so their velocity dispersions are likely underestimated. 


scaling relation. 
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TABLE 5 

Scaling Relations Between Dynamical Masses and SZ Signals 


Helation 

b 

a 

ay 

P{up\Y5ooD\) 
extended sample 
CIRS/HeCS 

P{Y500Dl\ap) 
extended sample 

U.l 1 D_q Q 22 
f) 1 Qi + 0.022 

U.iyi_ 0.022 

rv -| ■ 7 t:' + 0.029 

U.l / 0 _Q Q 3 Q 

z.OU_Q 29 

9 n 9 + 0.29 

z.uz 0.26 

O.uzu_o.oi9 

Q n 9 Q+ 0.021 

O.UZO_o.o 21 

q ni q+0.024 

O.UlO_o.o24 

"7 r-7+6.83 
'•O'_0.91 

—6 fi3+0-75 

0-Oo o.83 

u.uy/o_o.oo84 
nil 09 + 0 .0096 

u.iioz_0.0088 
on 4+0012 

n “^71 +0.035 
u.o/1_0.031 

U.oy4 0.026 

P{M20o\Msz) 

extended sample 

CIRS/HeCS 

P{^Sz\i^200) 

n ' 7 Q+U- 

<'5_o.12 

O'- <o_o 12 

'^-0.15 

1 '79+0.57 
-L-'^0.79 

14.0537";^^' 

14.006/Hn 

i4.uoy_Q Qgg 

- 24 . 2 / 1 I 3 ® 

u.z/y_Q 024 

0-308/^29 

0'-^o 0.i9 

P{(7p\Msz) 
extended sample 

P{MszWp) 

extended sample 

u.o±y_o 041 

0 33Q+0-043 
u.ooy_Q Q41 

1 03+^-12 
l.UO_Q Qg 

4-0«/_q q29 

Z.OOO_Q Q26 

_3 4-7+0.55 

-0.47 

9 q^+0.26 
z.ou o 34 

u.uy/o_Q QQgg 

n iiqq+0.0096 

U.liyO_Q QQgy 

U.ZUO_0.017 
n 999+0.017 
u.z^z Q 015 


Note. — Fits are of the relation P{y\x) assuming the linear form logy = a + blogx with intrinsic scatter criogy in the relation at fixed 
values of log x. 


4.2. Impact on the Tension Between Planck 
Cosmological Parameters from SZ versus CMB 

Cosmological constraints from Planck observations of 
the CMB predict a higher normalization of the clus¬ 
ter mass function (parameterized by flm and us) than 
the measured abundance from the SZ cluster detections 
(jPlanck Collaboration et al.ll2014all . As discussed in that 
paper, the cluster constraints are based on a scaling re¬ 
lation between SZ integrated Compton decrement and 
X-ray masses (calculated with the assumption of hydro¬ 
static equilibrium). They assume a hydrostatic mass 
bias due to non-thermal pressure support parameter¬ 
ized as Mtrue = (1 - b)MHSE where Mtme and Mhse 
are respectively the true cluster mass and the mass es¬ 
timated under the assumption of hydrostatic equilib- 
rium. Hydrodynamic simulations o f intracluster gas 
(jNagai et al.ll2007l : iNelson et al.ll20lJl predict a value of 
(1 — &) = 0.8, and the cosmological constraints are de¬ 
rived by allowing this parameter to vary in the range 
0 < 6 < 0.3. The tension between the SZ and CMB con¬ 
straints can be eliminated by assuming that the hydro- 
static mass bias is significant l y large r , (1—6) = 0.58±0.04 
(jPlanck Collaboration et ^ 12014a . 1201 ball . Note that 
the parameter b can have non-zero values either because 
of non-thermal pressure support or because of other cali¬ 
bration offsets (e.g., XMM-Newton temperature calib ra- 
tion, see [Israel et al.l[2M5t iSchellenberger et al.ll^OlSH . 

Estimates of hydrostatic mass bias from comparisons 
of X-ray and lensing mass estimates find smaller off¬ 
sets (e.g.. IVikhlinin et al.l l2009al: iMahdavi et al.l 120131 
lApplegate et al.ll2014ll . Recent revisions to systematic 
uncertainties in lensing mass estimates yield consistency 
in mass estim ates of individual clusters between different 
inves tigators (jvon der Linden et al.ll20l4lHoekstra et al.l 
nnn), yielding estimates of(l — &)~0.7 — 0.8, interme¬ 
diate between no hydrostatic mass bias and the large 
bias required to match CMB constraints. Alternatively, 
a new method of measuring weak lensing of the CMB by 
clusters yields 1/(1 — b) = 0.99 ± 0T9. cons istent with 
little to no mass bias (jMelin fc Bartlettll2014[ l. 

It is thus very interesting to see whether our dynamical 
mass estimates imply small hydrostatic mass bias (leav¬ 


ing tension between clusters and the CMB) or large hy¬ 
drostatic mass bias (alleviating tension between clusters 
and the CMB but aggravating tension among different 
cluster mass estimators). As mentioned in §3, our best- 
fit scaling relation is consistent with the Planck scaling 
relation based on hydrostatic mass estimates from XMM- 
Newton observations. Furthermore, Figure |9] shows that 
renormalizing this relation by assuming a hydrostatic 
mass bias of (1 — 6) = 0.58 ± 0.04 overpredicts the veloc¬ 
ity dispersion at fixed Ysz by an amount larger than the 
statistical uncertainties. That is, the CMB cosmologi¬ 
cal parameters predict significantly larger cluster velocity 
dispersions than our measured values. 

We next consider three possible explanations of the 
tension between the CMB normalization of the cluster 
mass scale and our measurement of the relation between 
velocity dispersion and integrated SZ decrement. First, 
we investigate whether X-ray selection (used for part of 
the sample at larger redshift) significantly impacts the 
resulting scaling relation. Second, we discuss the possi¬ 
bility of velocity bias (galaxies moving faster or slower 
than dark matter particles). Third, we discuss the pos¬ 
sible impact of massive neutrinos producing a smaller 
cluster abundance for a fixed matter power spectrum. 

4.3. Impact of Cluster Selection 

Scaling relations can be sensitive to the method of 
sample selection. The relation between dynamical mass 
and SZ signal could depend on whether the cluster sam¬ 
ple is selected from an optical catalog, an X-ray cata¬ 
log, or a SZ catalog, and whether the samples are flux- 
limited (detection-limited) or volume-limited. For in¬ 
stance, clusters with luminous cooling cores could be 
overrepresented in a flux-limited X-ray catalog com¬ 
pared to a mass-limited sample. The HeCS and CIRS 
cluster samples were drawn from X-ray-selected sam¬ 
ples. Thus, all clusters from these samples have mod¬ 
erately large X-ray fluxes. The Planck early release 
clusters contained several that were not previously de¬ 
tected in X-rays. Followup XMM-Newton observations 
of these clusters showed that they are in younger dy- 
namical states than the rest of the early release clusters 
((Planck Collaboration et ^12013^ . Thus, it is conceiv- 
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able that the HeCS and CIRS clusters are not a repre¬ 
sentative sample of Planck clusters. 

We use our Planck-selected sample to test whether 
the scaling relations depend on the selection technique. 
Specifically, we fit the scaling relations based on only 
clusters from the CIRS and HeCS samples, both of which 
are selected by X-ray flux. There is no significant change 
in the best-fit parameters for the X-ray selected sample 
compared to the SZ-selected sample (Table 00). Thus, 
the impact of X-ray selection versus SZ selection appears 
to be small, at least for the large and complete samples 
that we consider here. 

4.4. Can Velocity Bias Resolve the Tension? 

In numerical simulations, the velocity dispersion of 
randomly selected dark ma tter particles closely traces the 
mass of dark matter halos (lEvrard et al.l[2?)f)M . Observa- 
tionally, one measures the velocity dispersion of galaxies, 
which may move faster or slower than the underlying 
dark matter distribution. This “velocity bias” can be 
parametrized as hy = Ogxyj<J dm^ where <Jgxy and odm 
are the velocity dispersions of galaxies and dark matter 
particles respectively. 

If one assumes that the Planek CMB cosmological pa¬ 
rameters are correct, then the offset between the scal¬ 
ing relation we observe and the relation predicted by 
the CMB-based parameters provides information about 
the relation between galaxy dynamics and true cluster 
mass. In particular, significant negative velocity bias 
{by ~ 0.77) is required to bring the results into agree¬ 
ment. 

Modeling velocity bias in simulations is a very chal¬ 
lenging problem, due to both the uncertain physics in 
galaxy formation and evolution and the large dynamic 
range required to simulate individual cluster galaxies in 
a cosmological simulation. Some simulations follow the 
evolution of dark matter subhalos, but the evolution of 
galaxies may differ significantly because galaxies are ex¬ 
pected to form at the centers of dark matter subhalos 
and thus survive even after their dark matter halos are 
tidally stripped. 

Earlier, we used the consistency of the virial mass func¬ 
tion of X-ray selected clusters with cosmological con¬ 
straints from WMAP5, supernovae, and baryon acous¬ 
tic oscillations to conclu de that velocity b ias is small: 
cTgxy = (I.05±0.05)cr_DM (|Rines et al.ll2008fl . At present, 
there is no general agreement on the amount or even the 
sign of velocity bias, but the large negative velocity bias 
required for consistency with Planek CMB-based param¬ 
eters is not predicted by any current models. Simulations 
suggest that small samples of cluster galaxies restricted 
only to the brightest member s could be subject to neg¬ 
ative velocity bias of ^15% (|01d et al.ll20I3l : iWu et all 
l20Ri[l , but these simulations also suggest that large sam¬ 
ples such as the ones we analyze here should not be sub¬ 
ject to significant velocity bias. 

Note that a recent analysis of the redshift-space cor¬ 
relation function of high-mass galaxies from the SDSS 
Baryon Oscillation Sp ectroscopic Surv ey (BOSS) sug¬ 
gests that by Ks 0.86 (|Cuo et al.ll20I5[l . This negative 

^ Because Msz and YszD\ are closely related, we do not in- 
elude a separate fit for the o-p — Msz relation for the CIRS/HeCS 
subsample. 


velocity bias probably reflects the fact that even mas¬ 
sive and rich clusters contain very few high-mass galax- 
ies (e.g.. Figure 4 of I Guo et al.ll20I^ : thus, the mea¬ 
sured velocity bias is consistent with simulations that 
pre dict negative velocity bias for the brightest few galax- 


ies (lOld et al.l 120131: IWu et a 

.11201,811. It is also possible 

that the analysis of Guo et al. 

(I20I5I1 does not adequately 


mod el the impact of coherent inf all among satellite galax¬ 
ies (iHikaee fc Yamamotol I20I5I1 . Again, the spectro¬ 
scopic samples considered here include many galaxies be¬ 
low the characteristic absolute magnitude M* and are 
thus expected to have smaller bias than more luminous 
samples. 

Several recent simulations suggest that cluster galaxies 
should be positively biased by 5-15% depending on the 
detai l s of galaxy modelin g (|Lau et al.l 120101: iSaro et aD 
1201,81: iMunari et al.l 1201, 8fl . Positiye yelocity bias would 
further aggrayate the tension between the yelocity dis¬ 
persions we measure and the large SZ-mass normaliza¬ 
tion required to match the CMB data. 

Many of the simulations predict that the yelocity bias 
depends weakly on halo mass, so a more complete de¬ 
scription of velocity bias may require a virial scaling rela¬ 
tion with arbitrary slope (fixed by requires that the slope 
of the cTgxy — M relation is identical to the cjdm — M 
relation). Figure [H] shows several of these relations 
compared to our data, and Figure fTOl shows the param¬ 
eters of some of these models compared to the uncer¬ 
tainties in our observed scaling relation. The HeCS-SZ 
da ta and Planck rn asses are consistent wi th the models 
of iLau et al.l (1201011 . while the models of IMunari et al.l 
(l20I3ll lie far outside the observed relation. Importantly, 
although there is no consensus on the exact amount of 
velocity bias, none of the recent estimates are consis¬ 
tent with the large velocity bias {by ^ 0.77) required to 
reconcile the Planck SZ mass function with the CMB. 
Indeed, the discrepancy b etween our observed scaling re¬ 
lation and the models of IMunari et al.l (|20I3ll is in the 
opposite direction of the discrepancy required to reduce 
the CMB-SZ tension. 

4.5. Massive Neutrinos as a Solution? 

Neutrinos with significant masses can suppress the for¬ 
mation of large-scale structure. Thus, massive neutri¬ 
nos provide one possible explanation of the observed 
deficit of SZ clusters compared to the predictions from 
the best-fit ACDM model to the Planck CMB data. 
In particular, joint fits to CMB and SZ data from 
Planck yield estimates of total neutrino masses Sm,y = 
(0.40 ± 0.2I)eV when allowi ng the hydrostatic mass bias 
to var y between 0 and 0.3 (iPlanck Collaboration et ^ 
l20I4a|l . Adding baryon acoustic oscillation (BAO) mea- 
surements yields an estimate of Hm ,, = ( 0.20 ± 0.09)eV 
( Planck Collaboration et ^ l20I4ari . IWvman et al.l 
( 20I4f) point out that massive neutrinos not only alleviate 
tension between Planck CMB results and cluster abun¬ 
dance measurements, they also alleviate tension between 
Planck CMB results an d local measurements of the Hub¬ 
ble constant. Similarly. iMacCrann et aH (j20I5f l find that 
a similar tension exists between Planck CMB results and 
cosmic shear measurements; this tension can be partially 
alleviated with the introduction of a sterile neutrino. 
However, note that an alternate analysis of the cluster 
mass function using X-ray luminosities and weak lensing 
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Fig. 14.— Scaling relation between projected velocity dispersion 
tTp and the SZ mass proxy Msz based on the integrated Compton 
parameter Y^qqD^. The thick solid line shows the best-fit rela¬ 
tion of PicTplM^z) with the intrinsic scatter shown as the green 
band. The other lines show several predic tions of velocity bias . The 
dash-dotted line shows the predic tion ofIM unari et al.l II2013I ). the 
dotted line shows the predi ction ofILa u et al.l Il201fl l. the blue solid 
line shows the prediction of Saro^t_^ j20y) , the long-dashed line 
shows the prediction from Old'er^Tr||20T^ . and the blue short- 
dashed line show s the velocity bias for high-mass galaxies from 
I Guo et al.l II2015I1 . 


mass calibration yields reasonable agreement with the 
Planck C MB results, thus im plying no need for massive 
neutrinos (|Mantz et al.ll2015ll . The good agreement be¬ 
tween our measured velocity dispersions and those pre¬ 
dicted by the Planck SZ masses (assuming little velocity 
bias) supports the mass calibration used in the Planck SZ 
analysis. Our results therefore support the possibility of 
massive neutrinos as a solution to the CMB-SZ tension. 

Experimental measurements of neutrino oscillations 
place a lowe r limit of Hrr i y > 0.06eV [95% confi¬ 
dence level] (jCanozzi et al.l 1201411 . Constraints from 
the power spectrum of the Lyman-a forest from BOSS 
observations yield upper limits of Tirup <0.98 eV, 
or <0.16 eV when combined wit h Planck CMB data 
(|Palanaue-Delabrouille et al.ll2015[l . Thus, massive neu¬ 
trinos remain a plausible solution to the CMB-SZ ten¬ 
sion, but the required masses may produce tension with 
constraints from the Lyman-a forest. 

5. CONCLUSIONS 

The Planck satellite has produced a dramatic increase 
in the number of galaxy clusters with SZ mass estimates. 
Because the catalog includes many nearby clusters and 
covers the entire sky, many clusters in the Planck cata¬ 
log have existing mass estimates from galaxy dynamics. 
Here we measure 7,721 new redshifts in 24 clusters to 
obtain a large SZ-selected sample of 123 clusters with 
both dynamical and SZ mass estimates. To date, this 


is the largest sample of clusters used to compare veloc¬ 
ity dispersions and SZ mass estimates. We focus on a 
SZ-complete sample of 83 clusters. 

The measured velocity dispersions agree well with the 
predicted velocity dispersions from the cluster masses 
in the Planck SZ catalog and the virial scaling relation 
of dark matter particles. The cosmological parameters 
based on Planck CMB observations are not consistent 
with the mass function based on masses from the Planck 
SZ catalog. One way to resolve this tension is to al¬ 
low for mass bias in the SZ masses; large mass bias 
{Msz ~ 0.58Mtrue) is required to reconcile the CMB 
and SZ results. Such large mass bias is strongly disfa¬ 
vored by our results. 

In principle, velocity bias could allow galaxy veloc¬ 
ity dispersions to agree with the virial scaling relation 
for dark matter particles based on strongly biased SZ 
masses. However, no recent estimates of the amount of 
velocity bias are consistent with the large velocity bias 
{by « 0.77) required for this scenario. In fact, some mod¬ 
els of velocity bias have by > 1, a possibility that would 
further aggravate the tension between a possible SZ mass 
bias and our measured velocity dispersions. 

Departures from a standard ACDM cosmological 
model could resolve the tension between CMB and SZ 
cosmological parameter estimates. For example, signif¬ 
icant neutrino masses would decrease the amplitude of 
the power spectrum on clust er scales relative to the nor¬ 
maliza tion from the CMB (jPlanck Collaboration et Ml 
l2014a|l . In this scenario, Planck cluster masses could 
have little bias, and the excellent agreement between the 
measured velocity dispersions and the virial scaling rela¬ 
tion of dark matter particles would require that galaxy 
velocity bias is small (i.e., by ss 1). 

Future work on the equilibrium dynamics of cluster 
galaxies can test the possibility of large velocity bias: if 
large velocity bias is present, a Jeans analysis should re¬ 
veal that the cluster masses are larger than inferred by 
virial scaling relations (or by the caustic technique). Fu¬ 
ture simulations of the evolution of galaxies within clus¬ 
ters could test whether large velocity bias is plausible. If 
not, our results suggest that the tension between cosmo¬ 
logical parameters derived from CMB and SZ data may 
require extensions to the standard ACDM cosmological 
model. Observations of SZ-selected clusters at higher 
redshift could measure the evolution of cluster scaling 
relations and provide further insight into the origin of 
the CMB-SZ tension. 
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APPENDIX 

BAYESIAN PARAMETER ESTIMATION 


Define the likelihood p{D\0,M) the probability of measuring the set of data D when the model M is described by 
the set of parameters 6; the prior p{6\M) is the probability that the set 0 occurs. We are interested in estimating the 
probability density function (PDF) of the parameters 0 given our data set D 


p{e\D,M) 


p{D\e,M)p{e\M) 

p{D\M) 


(Al) 


Given the model M, we need to assume the likelihood p{D\9,M) and the prior p{9\M), whereas p{D\M) is a trivial 
normalization factor. 

In this work, we are interested in describing our data with linear correlations between pairs (A, Y) of the logarithm 
of the observables. In general, a number of unknown hidden variables produces a scatter in the linear correlation 
Y = a + bX. We model this scatter with a single parameter, the intrinsic dispersion (Tint. Therefore, given a 
measure Xi with uncertainty uxi, the probability of measuring with uncertainty cry. is p(Yi,aYi\9,Xi,aXi), where 
9 = {a, b, (Tint}- We assume the Gaussian likelihood 


p{D\9,M) 


n 


1 


exp 


' -{Y,-a-bX,f - 

2(t| 


(A2) 


where 

= O-fnt + ^Yi + • (A3) 

We assume independent flat priors for both a and b. For the intrinsic dispersion (Tint, which is positive defined, we 
assume 

p((Tint|M) =exp(-/ia;) (A4) 

r(r) 

where x = l/afnt; r(r) is the usual gamma function. This PDF describes a variate with mean r/p, and variance 
r /We set r = p = I0“® which guarantees an almost flat prior. 

To estimate the parameter PDF p{9\D, M), we perform a Markov Chain Monte Carlo (MCMC) sampling with 
the code APEMoST de veloped by Johannes Buchner and Michael Gruberbauer (jBuchner fc Gruberbauerri2nill : 
iGruberbauer et ^12009( 1. We obtain a fairly complete sampling with 2 x 10® MCMC iterations. The boundaries 
of the parameter space were set to [—100,100] for a and 6, and [0.01,100] for (Tint. The initial seed of the random 
number generator was set with the bash command GSL_RAND0M_SEED=$RANDDM. 

As the three best-fit parameters a, b, and (Tint of the Bayesian analysis, we adopt the medians derived from the 
posterior PDF p{9\D,M). Likewise, we adopt the boundaries of the 68% confidence levels around the medians as the 
uncertainties on these best-fit parameters. 


REFERENCES 


Adelman-McCarthy, J. K., et al. 2008, ApJS, 175, 297 
Ahn, C. P., et al. 2014, ApJS, 211, 17 
Andersson, K., et al. 2011, ApJ, 738, 48 
Andreon, S. 2010, MNRAS, 407, 263 
Applegate, D. E., et al. 2014, MNRAS, 439, 48 
Benson, B. A., et al. 2013, ApJ, 763, 147 
Biviano, A., & Girardi, M. 2003, ApJ, 585, 205 
Biviano, A., Murante, G., Borgani, S., Diaferio, A., Dolag, K., & 
Girardi, M. 2006, A&A, 456, 23 
Bocquet, S., et al. 2015, ApJ, 799, 214 
Bohringer, H., Burwitz, V., Zhang, Y.-Y., Schuecker, P., & 
Nowak, N. 2005, ApJ, 633, 148 
Bohringer, H., et al. 2000, ApJS, 129, 435 
—. 2004, A&A, 425, 367 

Bonamente, M., Joy, M., LaRoque, S. J., Carlstrom, J. E., Nagai, 
D., Sz Marrone, D. P. 2008, ApJ, 675, 106 
Bonamente, M., Landry, D., Maughan, B., Giles, P., Joy, M., Sz 
Nevalainen, J. 2013, MNRAS, 428, 2812 
Buchner, J., kc Gruberbauer, M. 2011, APEMoST (Automated 
Parameter Estimation and Model Selection Toolkit), 
http://apemost.sourceforge.net/ , commit from 2011-02-10 
Bullock, J. S., et al. 2UU1, MNRaS, (521, 559 

Capozzi, F., Fogli, G. L., Lisi, E., Marrone, A., Montanino, D., Sz 
Palazzo, A. 2014, Phys. Rev. D, 89, 093018 
Cypriano, E. S., Lima Neto, G. B., Sodre, Jr., L., Kneib, J.-P., &z 
Campusano, L. E. 2005, ApJ, 630, 38 
Czakon, N. G., et al. 2015, ApJ, 806, 18 
Danese, L., de Zotti, G., Sz di Tullio, G. 1980, A&A, 82, 322 
Diaferio, A. 1999, MNRAS, 309, 610 
—. 2009, ArXiv e-prints, arXiv/0901.0868 
Diaferio, A., Geller, M. J., & Rines, K. J. 2005, ApJ, 628, L97 


Ebeling, H., Edge, A. C., Bohringer, H., Allen, S. W., Crawford, 
C. S., Fabian, A. C., Voges, W., ^ Huchra, J. P. 1998, 

MNRAS, 301, 881 

Evrard, A. E., et al. 2008, ApJ, 672, 122 

Fabricant, D., Cheimets, P., Caldwell, N., Sz Geary, J. 1998, 

PASP 110 79 

Fabricant, D., et al. 2005, PASP, 117, 1411 
Faltenbacher, A., Sz Diemand, J. 2006, MNRAS, 369, 1698 
Geller, M. J., Diaferio, A., Rines, K. J., Sc Serra, A. L. 2013, ApJ, 
764, 58 

Gifford, D., Sz Miller, C. J. 2013, ApJ, 768, L32 

Gruberbauer, M., Kallinger, T., Weiss, W. W., & Guenther, D. B. 

2009, A&A, 506, 1043 
Guo, H., et al. 2015, MNRAS, 446, 578 

Hasselfield, M., et al. 2013, J. Cosmology Astropart. Phys., 7, 8 
Henry, J. P., Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, 
A. 2009, ApJ, 691, 1307 

Hikage, C., Sz Yamamoto, K. 2015, ArXiv e-prints 
Hoekstra, H., Hartlap, J., Hilbert, S., Sz van Uitert, E. 2011, 
MNRAS, 412, 2095 

Hoekstra, H., Herbonnet, R., Muzzin, A., Babul, A., Mahdavi, A., 
Viola, M., Sz Cacciato, M. 2015, MNRAS, 449, 685 
Hoekstra, H., et al. 2001, ApJ, 548, L5 

Hwang, H. S., Geller, M. J., Diaferio, A., Rines, K. J., Sz Zahid, 
H. J. 2014, ApJ, 797, 106 

Israel, H., Schellenberger, G., Nevalainen, J., Massey, R., Sz 
Reiprich, T. H. 2015, MNRAS, 448, 814 
Kurtz, M. J., Sz Mink, D. J. 1998, PASP, 110, 934 
Lau, E. T., Nagai, D., Sz Kravtsov, A. V. 2010, ApJ, 708, 1419 
MacCrann, N., Zuntz, J., Bridle, S., Jain, B., & Becker, M. R. 
2015, MNRAS, 451, 2877 













HeCS-SZ: Hectospec Survey of SZ-Selected Clusters 


13 


Mahdavi, A., Hoekstra, H., Babul, A., Bildfell, C., Jeltema, T., & 
Henry, J. P. 2013, ApJ, 767, 116 

Mamon, G. A., Biviano, A., & Boue, G. 2013, MNRAS, 429, 3079 
Mantz, A., Allen, S. W., Ebeling, H., Rapetti, D., & 

Drlica-Wagner, A. 2010a, MNRAS, 406, 1773 
Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010b, 
MNRAS, 406, 1759 

Mantz, A. B., et al. 2015, MNRAS, 446, 2205 
Marrone, D. P., et al. 2009, ApJ, 701, L114 
Melin, J.-B., & Bartlett, J. G. 2014, ArXiv e-prints 
Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005, 
ApJ, 623, L63 

Munari, E., Biviano, A., Borgani, S., Murante, G., & Fabjan, D. 
2013, MNRAS, 430, 2638 

Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007, ApJ, 655, 98 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Nelson, K., Lau, E. T., & Nagai, D. 2014, ApJ, 792, 25 
Old, L., Gray, M. E., & Pearce, F. R. 2013, MNRAS, 434, 2606 
Old, L., et al. 2014, MNRAS, 441, 1513 

Palanque-Delabrouille, N., et al. 2015, J. Cosmology Astropart. 
Phys., 2, 45 

Perrott, Y. C., et al. 2014, ArXiv e-prints 
Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & 
Melin, J.-B. 2011, A&A, 534, A109 
Planck Collaboration et al. 2011, A&A, 536, All 
—. 2013a, A&A, 550, A130 
—. 2013b, A&A, 550, A131 
—. 2014a, A&A, 571, A20 
—. 2014b, A&A, 571, A29 


—. 2015a, ArXiv e-prints 

—. 2015b, ArXiv e-prints 

Rines, K., & Diaferio, A. 2006, AJ, 132, 1275 

Rines, K., Diaferio, A., & Natarajan, P. 2007, ApJ, 657, 183 

—. 2008, ApJ, 679, LI 

Rines, K., Geller, M. J., & Diaferio, A. 2010, ApJ, 715, L180 
Rines, K., Geller, M. J., Diaferio, A., & Kurtz, M. J. 2013, ApJ, 
767, 15 

Rines, K., Geller, M. J., Kurtz, M. J., & Diaferio, A. 2003, AJ, 
126, 2152 

Rozo, E., et al. 2010, ApJ, 708, 645 
Ruel, J., et al. 2014, ApJ, 792, 45 

Saro, A., Mohr, J. J., Bazin, G., & Dolag, K. 2013, ApJ, 772, 47 
Schellentierger, G., Reiprich, T. H., Lovisari, L., Nevalainen, J., & 
David, L. 2015, A&A, 575, A30 
Serra, A. L., & Diaferio, A. 2013, ApJ, 768, 116 
Serra, A. L., Diaferio, A., Murante, G., & Borgani, S. 2011, 
MNRAS, 412, 800 
Sifon, C., et al. 2013, ApJ, 772, 25 

Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comments on 
Astrophysics and Space Physics, 4, 173 
Vikhlinin, A., et al. 2009a, ApJ, 692, 1033 
—. 2009b, ApJ, 692, 1060 

von der Linden, A., et al. 2014, MNRAS, 443, 1973 

Wu, H.-Y., Hahn, O., Evrard, A. E., Wechsler, R. H., & Dolag, 

K. 2013, MNRAS, 436, 460 

Wyman, M., Rudd, D. H., Vanderveld, R. A., & Hu, W. 2014, 
Physical Review Letters, 112, 051302 
Zwicky, F. 1933, Helv. Phys. Acta, 6, 110 
—. 1937, ApJ, 86, 217 



